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Abstract —In this paper, we introduce a sophisticated path loss 
model Incorporating hoth llne-of-sight (LoS) and non-line-of-slght 
(NLoS) transmissions to study their impact on the performance of 
dense small cell networks (SCNs). Analytical results are obtained 
for the coverage prohahility and the area spectral efficiency 
(ASE), assuming hoth a general path loss model and a special 
case with a linear LoS prohahility function. The performance 
impact of LoS and NLoS transmissions in dense SCNs in terms 
of the coverage prohahility and the ASE is significant, hoth 
quantitatively and qualitatively, compared with the previous work 
that does not differentiate LoS and NLoS transmissions. Our 
analysis demonstrates that the network coverage prohahility first 
increases with the increase of the hase station (BS) density, and 
then decreases as the SCN becomes denser. This decrease further 
makes the ASE suffer from a slow growth or even a decrease 
with network densification. The ASE will grow almost linearly as 
the BS density goes ultra dense. For practical regime of the BS 
density, the performance results derived from our analysis are 
distinctively different from previous results, and thus shed new 
insights on the design and deployment of future dense SCNs. 

Index Terms —stochastic geometry. Homogeneous Poisson Point 
Process (HPPP), Line-of-Sight (LoS), Non-Llne-of-Sight (NLoS), 
dense small cell networks (SCNs), coverage probability, area 
spectral efficiency (ASE). 

I. Introduction 

Driven by a new generation of wireless user equipment 
(UE) and the proliferation of bandwidth-intensive applications, 
mobile data traffic and network load are increasing in an 
exponential manner, and are straining current cellular networks 
to a breaking point [1]. In this context, small cell networks 
(SCNs), comprised of remote radio heads, metrocells, pico- 
cells, and/or femtocells, have attracted much attention as one 
of the most promising approaches to increase network capacity 
and meet the ever-increasing capacity demands [2], [3]. SCNs 
can significantly enhance network capacity through a high 
spatial spectrum reuse, e.g., network capacity could potentially 
grow linearly with the number of small cells [2]. Due to such 
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capacity gains, the orthogonal deployment of SCNs within the 
existing macrocell network, i.e., small cells and macrocells 
operate on different frequency spectrum (Small Cell Scenario 
#2a defined in [3]), have gained much momentum in the 
design of the 4th-generation (4G) Long Term Evolution (LTE) 
networks by the 3rd Generation Partnership Project (3GPP). 
Orthogonal deployments of dense SCNs are also envisaged as 
the workhorse for capacity enhancement in the 5th-generation 
(5G) networks, aided by its easy deployment arising from 
its low interaction with the macrocell tier, e.g., no inter-tier 
interference [2]. In this paper, our focus is on these orthogonal 
deployments of dense SCNs. 

In order to deploy dense SCNs in a cost-effective manner, 
vendors and operators need foremost a deep theoretical under¬ 
standing of the implications that small cells bring about. Being 
aware of the need for such knowledge, the wireless industry 
and research community have been working relentlessly on 
the modeling and the analysis of the dense SCN deployments. 
However, up to now, most studies on SCNs have considered 
only simplistic path loss models that do not differentiate Line- 
of-Sight (LoS) and Non-Line-of-Sight (NLoS) transmissions 
[4-7]. It is well known that LoS transmission may occur when 
the distance between a transmitter and a receiver is small, 
and NLoS transmission is common in office environments 
and in central business districts. Moreover, when the distance 
between a transmitter and a receiver decreases, the probability 
that a LoS path exists between them increases, thereby causing 
a transition from NLoS transmission to LoS transmission 
with a higher probability. In this paper, we will study the 
performance impact of such NLoS-to-LoS transition in dense 
SCNs. The main contributions of this paper are as follows: 

• Analytical results are obtained for the network perfor¬ 
mance of the coverage probability and the area spectral 
efficiency (ASE), using a proposed general path loss 
model incorporating both LoS and NLoS transmissions. 

• Using the above results, numerically tractable integral- 
form expressions for the coverage probability and the 
ASE are further obtained for a 3GPP path loss model 
with a linear LoS probability function. Our analysis can 
be readily extended to deal with more complicated path 
loss models by approximating the corresponding LoS 
probability function as a piece-wise linear function. 

• Our theoretical analysis reveals an important finding, i.e., 
when the density of small cell base station (BSs) is larger 
than a threshold, the network coverage probability will 
decrease as small cells become denser, which in turn 
can make the ASE suffer from a slow growth or even a 
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decrease. Thereafter, the ASE will grow almost linearly 
as the small cell BS density increases above another 
threshold, larger than the previous one. This finding that 
the ASE performance may suffer from a slow growth 
or even a decrease as the BS density increases, is not 
only quantitatively but also qualitatively different from 
previous study results with a simplistic path loss model 
that does not differentiate LoS and NLoS transmissions. 
Although our conclusion is made from the investigated 
set of parameters recommended by the 3GPP for SCNs, 
which sheds valuable insights on the design and deploy¬ 
ment of practical SCNs, it is of significant interest to 
further study the generality of our conclusion in other 
network models and with other parameter sets. 

The remainder of this paper is structured as follows. Sec¬ 
tion II provides a brief review on stochastic geometry and 
compares the closest related works to our work. Section III 
describes the system model. Section IV presents our main 
analytical results on the coverage probability and the ASE, fol¬ 
lowed by their application in a 3GPP study case addressed in 
Section V. The numerical results are discussed in Section VI, 
with remarks shedding some new light on the performance 
and the deployment of dense SCNs. Einally, the conclusions 
are drawn in Section VII. 

II. Related Work 

In stochastic geometry, BS positions are typically modeled 
as a Homogeneous Poisson Point Process (HPPP) on the plane, 
and closed-form expressions of coverage probability can be 
found for some scenarios in single-tier cellular networks [4] 
and multi-tier cellular networks [5], [6]. A general treatment of 
stochastic geometry can be found in [7]. The major conclusion 
in [4-7] is that neither the number of cells nor the number 
of cell tiers changes the coverage probability in interference- 
limited fully-loaded wireless networks. However, these works 
consider a simplistic path loss model that does not differentiate 
LoS and NLoS transmissions. In contrast, in this paper, we 
consider a more complete path loss model incorporating both 
LoS and NLoS transmissions to study their impact on the 
performance of dense SCNs. 

Notions that are similar to LoS and NLoS transmissions 
have been previously explored in the building blockage study 
in [8] and the indoor communication network in [9]. In [8], 
the authors proposed a microscopic performance analysis 
framework to model the random blockage effect of buildings, 
and analyze its impact on cellular network performance. 
Lurther refinement and verification of the proposed model 
in [8] is needed, especially to consider reflections which are 
an important contributor to coverage in urban areas. In [9], 
the authors present an analytical study of indoor propagation 
through walls, and showed that the throughput does not scale 
linearly with the density of small cells. Different from [9], in 
this paper, we investigate outdoor dense SCNs. 

The closest related works to the one in this paper 
are [10], [11] and [12]. 

In [10], the authors assumed a multi-slope piece-wise path 
loss function. Specifically, assuming that the distance between 


a BS and a UE is denoted by r in km, then the path loss 
associated with distance r can be formulated as 


c (r) = < 


Cl {r) , 

C2 (t) , 


when 0 < r < di 
when c?i < r < c ?2 


, Cn (r), when r > 


( 1 ) 


where the path loss function ( (r) is segmented into N pieces, 
with each piece and each segment break point denoted by 
Cn (r) and dn, respectively. 

In [11], the authors treated the event of LoS or NLoS 
transmission as a probabilistic event for a millimeter wave 
communication scenario. Specifically, the path loss associated 
with distance r is formulated as 


cw 


(r) , with probability Pr^ (r) 

(r) , with probability (l — Pr^ (r)) 


where (r), C^^ (t) and Pr^ (r) are the path loss function for 
the case of LoS transmission, the path loss function for the 
case of NLoS transmission and the LoS probability function, 
respectively. To simplify the analysis, the LoS probability 
function Pr^ (r) was approximated by a moment matched 
equivalent step function in [11]. 

In [12], the authors used the same path loss model as 
in (2) and considered the approximation of Pr^ (r) as an 
exponentially decreasing function. The results in [12] are less 
tractable than those in [10] and [11]. This is because the 
exponentially decreasing LoS probability function, albeit more 
practical than the step function in [11], is still difficult to deal 
with in the analysis. 

In this paper, we extend the works in [10-12] to an 
even more general scenario to improve the following aspects. 
In [10], the multi-slope piece-wise path loss model in (1) does 
not fit well with the practical model defined by the 3GPP, 
in which the path loss function is not a one-to-one mapping 
to the distance. In [11], the single-piece path loss model 
and the proposed step function are not compatible with the 
practical piece-wise path loss functions, the detailed modeling 
of which is presented in Section III. In [12], the considered 
path loss model is very practical. However, the generality 
of the study and the tractability of the analysis need to be 
improved. Compared with [10-12], the novelties of our paper 
are summarized as follows. 


• Extending the existing works to a more general 
scenario: We propose a path loss model that features 
piece-wise path loss functions with probabilistic LoS and 
NLoS transmissions. The proposed path loss model, to 
be presented in Section III, is general and can be applied 
to several channel models that capture LoS and NLoS 
transmissions [10-14]. Moreover, new analytical results 
are obtained for the coverage probability and the ASE. 

. Presenting more tractable results for a special path 
loss model with a linear LoS probability function: We 
derive numerically tractable integral-form expressions for 
the coverage probability and the ASE for a 3GPP path 
loss model with a linear LoS probability function. The 
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results are more simple and tractable than those in [12], 
and no approximation such as the step function [11] is 
used in our analysis. Although the results in [10] are 
superior than ours in terms of tractability, our analysis 
can characterize practical networks more accurately due 
to the consideration of probabilistic LoS events and a 
linear LoS probability function. Note that the inclusion of 
a linear LoS probability function is not trivial compared 
with [10] in the sense that, in our model a UE may be 
associated with a BS that is not the nearest BS to the UE, 
but where such BS has a LoS path to the UE resulting 
in the smallest path loss (i.e., with the largest ( (r)). In 
contrast, in [10], a UE always connects with its nearest 
BS. Moreover, our study on the path loss model with 
the linear LoS probability function can be extended to 
various cases. In this paper, we show that the approach 
of approximating a LoS probability function by a piece- 
wise linear function and invoking our results, can deal 
with complicated path loss models in a tractable manner. 

• Disclosing a new finding on the ASE performance: We 
present a new finding that the ASE performance without 
the assumption of near-field path loss exponents may 
suffer from a slow growth or even a decrease with BS 
densification. This new finding is not only quantitatively 
but also qualitatively different from the results given 
by [10-12]. Note that in [10], a deterministic rate based 
on the signal to interference plus noise ratio (SINK) 
threshold is assumed for the typical UE, no matter what 
the actual SINK value is. Our definition of the ASE, to 
be formally presented in Section III, considers a more 
realistic SINR-dependent rate, which leads to a more 
complex analysis that requires one more fold of numerical 
integral compared with [10]. 


III. System Model 


We consider a downlink (DL) cellular network with BSs 
deployed in a plane according to an HPPP <I> of intensity 
A BSs/km^. UEs are Poisson distributed in the considered 
network with an intensity of A™ UEs/km^. Note that A^^ 
is assumed to be sufficiently larger than A so that each BS 
has at least one associated UE in its coverage. As in (1) and 
(2), the distance between an arbitrary BS and an arbitrary UE 
is denoted by r in km. Considering practical LoS and NLoS 
transmissions, we propose to model the path loss with respect 
to distance r as (3), which is shown on the top of next page. 

In (3), the path loss function ({r) is segmented into N 
pieces with each piece denoted by Cn (t). Besides, (r), 
(r) and Pr^ (r) are the n-th piece of path loss function for 
the LoS transmission, the n-th piece of path loss function for 
the NLoS transmission, and the n-th piece of the LoS prob¬ 
ability function, respectively. Note that the proposed model 
is general and can be applied to several channel models that 
capture LoS and NLoS transmissions [10-14]. 

Moreover, (r) and (r) in (3) are modeled as 


Cn (r) 


=ALr-“n, for LoS 
(Nl ^ ^ for NLoS ’ 


( 4 ) 


where and A^^,n € {1, 2,..., W} are the path losses at 
a reference distance r = 1 for the LoS and the NLoS cases 
in Cn (r), respectively, and and a^^,n G {1,2,...,A^} 
are the path loss exponents for the LoS and the NLoS cases 
in („ (r), respectively. In practice, Aj^, and are 

constants obtained from field tests [13], [14]. Eor convenience, 
{CnW} and {Cr(r)} are further stacked into piece-wise 
functions written as 


c 


Path 


(r) = < 


/•Path 

51 

/•Path 

52 


h) , 


when 0 < r < di 
when di < r < d 2 


( 5 ) 


when r >(iAr_i 


where the string variable Path takes the value of “L” and 
“NL” for the LoS and the NLoS cases, respectively. 

In (3), Pr^ (r) ,n G {1,2,..., N} is the n-th piece proba¬ 
bility function that a transmitter and a receiver separated by a 
distance r has a LoS path, which is typically a monotonically 
decreasing function with r. Eor convenience, {Pr^(r)} is 
further stacked into a piece-wise LoS probability function 
expressed as 


Pr\ (r), 

when 0 < r < di 

PrL(r), 

when di < r < d 2 

P4 (t) , 

when r > d^v-i 


Our model is consistent with the ones adopted in the 
3GPP [13], [14]. It should be noted that the considered path 
loss model shown in (3) includes the models in [10-12] as its 
special cases. More specifically, for [10], all the LoS probabil¬ 
ities are zero, i.e., assuming Pr^ (r) = 0,VnG{l,2,...,iV} 
in (3). Eor [11] and [12], there is only one piece of path loss 
function with one LoS path loss exponent and one NLoS path 
loss exponent, i.e., assuming = 1 in (3). 

Another important note is on the practical usage of the pro¬ 
posed path loss model, which is summarized in the following: 

« As have been addressed in [10], the standard single¬ 
slope path loss function does not accurately capture 
the dependence of the path loss exponent on the link 
distance in many future networks, such as the multi¬ 
ray transmission environments, the dense/clustered small 
cells, the millimeter wave communications with blockage, 
etc. Therefore, a dual or more slope path loss function 
such as (5) should be considered. 

« The realistic LoS probability functions usually take com¬ 
plicated mathematical forms, e.g., in the 3GPP stan¬ 
dards [13], which will be addressed in more detail in 
Section V. Therefore, to achieve both analytical tractabil¬ 
ity and result accuracy, it is desirable to approximate the 
complicated LoS probability function as a few pieces of 
elementary functions, e.g., linear functions. Such piece- 
wise LoS probability function is well captured by (6) and 
can be studied in our framework. 

• With the justification of both (5) and (6), we can conclude 
that the proposed path loss model in (3) is versatile to 
cope with various practical network scenarios. 
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Cl {r) 


C(r) = < 


C 2 (r) 


Cl (r), with probability Pri (r) 

C 5 ^^ (r), with probability (l—Pri(r))’ 
C 2 (’’) I with probability Pr^ (r) 

C^^ (r), with probability (l —Pr 2 (r'))’ 


when 0 < r < di 
when di < r < 6,2 


Cn (r) 


Cjc (^): 
crw, 


with probability Pr^ (r) 
with probability (l — Pr^ (r))’ 


when r > dN-i 


( 3 ) 


In this paper, we assume the following user association 
strategy (UAS). Each UE should be associated with the BS 
having the smallest path loss (i.e., with the largest C {r)) to 
the UE [11], [12]. Note that in our previous work [15] and 
some existing works [4], [10], it was assumed that each UE 
should be connected to the BS with the closest proximity. Such 
assumption is not appropriate because in practice it is possible 
for a UE to associate with a BS that is not the closest one but 
the one with the minimum path loss. 

Einally, we assume that each BS/UE is equipped with an 
isotropic antenna, and as a common practice in the field [4- 
6,10,12], the multi-path fading between an arbitrary BS and an 
arbitrary UE is modeled as independently identical distributed 
(i.i.d.) Rayleigh fading. 


IV. Analysis for the Proposed Path Loss Model 


Using the property of the HPPP, we study the performance 
of SCNs by considering the performance of a typical UE 
located at the origin o. We first investigate the probability that 
the typical UE is covered by its associated BS. This coverage 
probability is defined as the probability that the UE’s SINR is 
above a per-designated threshold 7 : 

p“''(A, 7 ) = Pr[SINR> 7 ], (7) 


where the SINR is computed by 


PC (t) h 




Ir + Nq 


( 8 ) 


where h is the channel gain of Rayleigh fading, which is 
modeled as an exponential random variable (RV) with the 
mean of one, and P and Nq are the transmission power 
of each BS and the additive white Gaussian noise (AWGN) 
power at each UE, respectively. Moreover, N is the aggregate 
interference given by 


/, = ^ P^,g„ (9) 

i: hi^^\ho 


where ho is the BS serving the typical UE, which is located 
at distance r from the typical UE, and hi, jdi and gi are the 
j-th interfering BS, the path loss associated with hi and the 
multi-path fading channel gain associated with hi, respectively. 

Moreover, according to [12], the area spectral efficiency 
(ASE) in bps/Hz/km^ for a given A can be defined as 

pOO 

= X log2(l-f 7)/r(A,7)d7, (10) 

7 70 


where 70 is the minimum working SINR for the considered 
SCN, and fr (A, 7 ) is the probability density function (PDE) 
of SINR observed at the typical UE for a particular value of 
A. Note that the ASE defined in this paper is different from 
that in [ 10 ], where a deterministic rate based on 70 is assumed 
for the typical UE, no matter what the actual SINR value is. 
The ASE definition in (10) is more realistic due to the SINR- 
dependent rate, but it is more complex to analyze, as it requires 
one more fold of numerical integral compared with [ 10 ]. 


Based on the definition of p^°'^ (A, 7 ) in (7), which is the 
complementary cumulative distribution function (CCDE) of 
SINR, fr (A, 7 ) can be computed by 


fr (A, 7 ) 


g(l-p-°MA, 7 )) 

97 


( 11 ) 


Given the definition of the coverage probability and the ASE 
respectively presented in (7) and (10), in the following we 
will analyze the two performance measures for the considered 
UAS. Based on the proposed path loss model in (3), we present 
our main result on p^°'' (A, 7 ) in Theorem 1. 


Theorem 1. Considering the path loss model in (3), p™'’ (A, 7 ) 
can be derived as 


N 




( 12 ) 


Pr 


PCir)h ; 
Ir+No P I 


n—1 

Ir + No 


> 7 


fkn(.r)dr, = 


f^n (^) do and d^ are re- 


where Tf = 

spectively defined as 0 and 00 . Moreover, ^ (r) and fr\ (r) 

are given by 

fkn (.r) = exp 1 ^- ^ (1 - M (u)) 2TruXdv^ 

X exp ( — / Pr^ (u) 2'KuXdu ) 


xPr^ {f') X 27rrA, (d„_i < r < d„), (13) 

and 

fR% (r) = exp J Pr^ {u) 2 TTuXdu^ 

X exp J (1 — Pr^ (u)) 2 ‘nuXd’u^ 

X (1 - (r)) X 27rrA, (d„_i < r < dn ), (14) 
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where ri and r 2 are determined by 

ri = arg{C'^^(ri) = (r)} , 

ri 

and 


r2 = arg{C^(r2) = W} ■ 


(15) 

(16) 


Furthermore, Pr 


_ Ir + No > ^ 

respectively computed by 


r2 

PCir)h 


and Pr 


EC(r)h ^ 
Ir+No ^ ' 


Pr 

and 

Pr 


PCk{r)h 

Ir + Nr) 


> 7 


=exp 


PCHr)h 

Ir + No 


> 1 


=exp 


jNq 

PCk (r) 

yNo 


^Ir 


PCk (r) 


, (17) 


^7 


7 


pe(o; \pc^Hr)j’ 

(18) 

where (s) is the Laplace transform of N evaluated at s. 


Proof: See Appendix A. 


As can be observed from Theorem 1, the piece-wise path 
loss function for LoS transmission {Cra(^)}> the piece-wise 
path loss function for NLoS transmission (^)}> the 
piece-wise LoS probability function {Pr^(r)} play active 
roles in determining the final result of p“''(A, 7 ). We will 
investigate their impacts on network performance in detail in 
the following sections. 

Plugging p“''(A, 7 ) obtained from (12) into (11), we can 
get the result of the ASE using (10). 

Regarding the computational process to obtain (A, 7 ) 
presented in Theorem 1, for a general case, three folds 
of integrals are respectively required for the calculation of 

and where the 

string variable Path takes the value of “L” (for the LoS case) 
or “NL” (for the NLoS case). Note that an additional fold of 
integral is needed in ( 10 ) for the calculation of (A, 70 ), 
making it a 4-fold integral computation. 


V. Study of a 3GPP Special Case 


As a special case for Theorem 1, we consider the following 
path loss function, ( (r), adopted by the 3GPP [13], i.e.. 


C(t) 


A^r , with probability Pr^ (r) 

with probability (l — Pr^ (r)) ’ ^ ^ 


together with a linear LoS probability function of Pr^ (r), also 
adopted by the 3GPP [14], i.e.. 


Pr^ (r) 


1 - 77 , 0 <r<di 
0 , r > di 


( 20 ) 


where di is a parameter that determines the decreasing slope 
of the linear function Pr^ (r). 

Considering the general path loss model presented in (3), 
the path loss model shown in (19) and (20) can be deemed as 
a special case of (3) with the following substitution: N = 2, 
Cl- (r) = (r) = ALr-“^ Cr (r) = (^) = 

Prl; (r) = 1 — and Pr^ (r) = 0. For clarity, this 3GPP 
special case is referred to as 3GPP Case 1 in the sequel. 


It should be noted that 3GPP Case 1 is compatible with 
dense SCNs, because both the exponential path loss function 
in [13] and the LoS probability function in [14] are valid for 
small cells only (referred to as microcells in [14]). It should 
also be noted that there is another LoS probability function 
defined in [13], which takes a more complicated mathematical 
form given by 

Pr^ (r) = 0.5 — min < 0.5,5 exp 




where Ri and R 2 are shape parameters to ensure the continuity 
of Pr‘-(r). 

To show how Pr^ (r) in (21) can be fitted into our proposed 
path loss model, we can re-formulate ( 21 ) as 


Pr^ (r) 


1 — 5 exp {—Ri/r ), 

5 exp(-r/i? 2 ), 


where di 


Ri 

Into' 


0 < r < di 
r > di 


( 22 ) 


The combination of the LoS probability function in (22) and 
the path loss function in (19) can then be deemed as a special 
case of the proposed path loss model in (3) with the following 
substitution: N = 2, (r) = C 2 (^) = (t) = 

(t) = Pr^ (r) = 1 — 5 exp (—i?i/r), and 

Pr 2 (r) = 5 exp (—r/i? 2 )- For clarity, this combined case with 
both the path loss function and the LoS probability function 
coming from [13] is referred to as 3GPP Case 2 in our paper. 
Note that 3GPP Case 2 was treated in [12] by approximating 
Pr^ ir) in (21) as an exponentially decreasing function. How¬ 
ever, the results in [ 12 ] are less tractable than ours as well as 
those in [ 10 ] and [ 11 ], because the approximated function is 
still difficult to deal with in the theoretical analysis. 


In the following, we first investigate 3GPP Case 1 in our 
case study, because it gives more tractable results than 3GPP 
Case 2 treated in [12], as will be shown in the following 
subsections. Thereafter, we will numerically investigate 3GPP 
Case 2 using Theorem 1 in Section VI, and will show that 
similar conclusions like those for 3GPP Case 1 can also be 
drawn for 3GPP Case 2. More importantly, we will extend 
3GPP Case 1 to study an approximated 3GPP Case 2 in 
Section VI, thus showing the usefulness of studying 3GPP 
Case 1 with the linear LoS probability function. 


To sum up, taking the linear LoS probability function 
from [14] to create 3GPP Case 1 not only allows us to obtain 
more tractable results, but also help us to deal with more 
complicated path loss models in practice. 

For 3GPP Case 1, according to Theorem 1, p“'' (A, 7 ) can 
then be computed by 

2 

P^°^X,l) = J2iTn+Tn^)- ( 23 ) 

n—1 


In the following subsections, we investigate the computation 
of T|-, T 2 , and respectively. 
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A. The Computation ofT^ 

From Theorem 1, Tf can be obtained as 




‘I (24) 

where Ci (^) = 24'^'r““ from (19) is plugged into the step 
(a) of (24) and (s) is the Laplace transform of RV Ir 
evaluated at s. 

In (24), according to Theorem 1 and (20), we have 

fk,i (^) = exp J \^2TTudu^ 


= exp I —TrXr^ + 2 ttX 


3di 3di 


X ( 1 — — ) 27rrA, (0 < r < di ), 


(25) 


where ri = 


i^y 


according to (15). 


Besides, to compute ^ j in (24) for the range of 
0 < r < c?i, we propose Lemma 2. 


Lemma 2. | -pjr ) hi the range of 0 < r < di can be 

calculated by 


■^L 




PA^ 


exp^-27rA|pi|Q;^, 1, 1, , 

X exp(^^pi||r^,2, (7r“'") , 2, 


27rA 


X exp ——Pi 2, I —n—r 


NL O 

I 


27rX 


H 


/ 7 ANL 


AL 


,di 

-1 


,n 


X exp|-27rAp2(a^^, 1, ) ,di)), (0 < r < di) 




where 

Pi (a,/3,f,d) = 
and 

P2 (a,/3,f,d) = 


d(/3+i) 


2 F 1 


/3 + 1 

d-(“-^-i) 




a 


t{a — /3 — 1) 


X 2d^l 




0 ) 


,{a>P+l), (28) 


a a td° 

where 2 F 1 is the hyper-geometric function [16], 


Proof: See Appendix B. 

To sum up, Tf can be evaluated as 




7^\ 

paA 


where /)j 1 (r) and j given by (25) and (26), 

respectively. 

B. The Computation of 

From Theorem 1, can be obtained as 


y NL ^ 


(a) 




-1 
yNo 




Jo (~ ^P7LN^° ) ( ?7Nt:) fpi 0 0 


PCr (r) 

C ^ nl \ 
p^NL 


fkki^)dr 


(30) 


where (r) = “ from (19) is plugged into the step 

(a) of (30). 

In (30), according to Theorem 1 and (20), we have 

/S (^) = ®xp J APr^ {u) 27rudit^ 

X exp r A (1 — Pr^ (u)) 2 TTudii^ 

X 27rrA, (0<r<di), (31) 

where r 2 = ( j according to (16). Since the nu¬ 

merical relationship between r 2 and di affects the calculation 
of the first multiplier in (31), i.e., exp (— 0 APr^ (u) 2-Kudu'), 
we should discuss the cases of 0 < r 2 < di and r 2 > di. 


have 


If 0 < r 2 < di, i.e., 0 < r < pi = d0 

fkk (r) = exp (^- J ^ 2 TTudu^ 

-2'Kudu\ ( — I 2'KrX 

1 J \dij 


we 


X exp I — / A 


(26) 


= exp|-,A.,+2,A( 3 ^- 5 ^ 


X I — ) 27rrA, (0 < r < pi). 


(32) 


Otherwise, if r 2 > di, i.e., yi < r < di, we have 

ffti ^ ~ 

X exp ( — / X—2'Kudu | ( ^ | 2'KrX 

\ Jo di J \dij 


= exp 


TrAdf 27rAr^ A 

3 3^ J 


X ( — ) 27rrA, {yi < r < di). (33) 
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Besides, to compute ^ j in (30) for the range 

of 0 < r < di, we propose Lemma 3 in the following. 
Note that since the calculation of (r) is divided into two 
cases respectively shown in (32) and (33), the calculation of 

/ cNL \ 

■^ir I 1 should also be divided into those two cases, be¬ 

cause the first case includes both LoS and NLoS interference, 
while the second case only considers NLoS interference. 

Lemma 3. .jSf/,, I 1 in the range of 0 < r < di can be 

divided into two cases, i.e., 0 < r < yi and yi < r < di. The 
results are as follows. 




"fr 


I p^NL 


exp —27rApi 1, 


7^^ 


, di 


-f27rApi la , 1, I 




TlNL 


,r 2 


( 2'k\ / T r, / 1-^ 


SttA ( T „ f 1 -^ 


di 

-1 


A 2 


X exp 


27rA 
di 
2'kX 


Pi ^a^^,2, ,di^ 

a-,2,(7r“^y\r)) 

X exp ^-27rAp2 ^a^^, 1, , di^ ^ , 


and 


/ \ 

/ IT \ 


Tr= / 

Jo 


exp -- 


NL , \ 

7r“ No\ 

p^NL 




/ NL 

(1 .NL 


ip^NL 


mr) 


0 < r < 2/1 


dr 


7 “ 


exp -- 






p^NL 

mr) 


/ NL’ 

ITT 1 .NL 


2/1 < r < di 


dr, (36) 


^p^NL 

where (32), (33), (34) and (35) are plugged into (36). 

C. The Computation ofT^ 

From Theorem 1, Tf can be derived as 

According to Theorem 1 and (20), (r) is given by 

/fl,2 (^) ^ (“ / ^ (1 “ (“)) 27rMdu^ 

X exp J APr^ (tt) 27ritdu^ x 0 x 27rrA 
= 0, (r>di). (38) 


Plugging (38) into (37), yields 

Tt = 0 . 


(39) 


D. The Computation of 

From Theorem 1, can be derived as 


(0 < r < 2/i), (34) 


I^pt^nl j 

exp 2, , di^ 

xexp ^-27rAp2 ^a^^,l, (^ 7 r“ ^ , 

( 2/1 < r < di), (35) 

where pi (a, jd, t, d) and p 2 (a, jd, t, d) are defined in (27) and 
(28), respectively. 

Proof: See Appendix C. ■ 

To sum up, can be evaluated as 

rVl 


(40. 

where (r) = from (19) is plugged into the step 

(a) of (40). 

In (40), according to Theorem 1 and (20), we have 
/fl !2 (?^) =exp (- J ^ 2 TTudu\ 


X exp ( — / \—2TTudu — / X2'Kudu ) 2'KrX 
\ Jo di Ja, J 

= exp (—7rAr^) 27rrA, (r > di). (41) 

Besides, to compute .if/,, I ) in (40) for the range of 

r > di, we propose Lemma 4. 


Lemma 4. .if/,, ^/f^TrrJ in the range of r > di can be 
calculated by 

(-2^Ap2 (^aNL,l, (7r“'''') , 

(r>di), (42) 

where p 2 (a,jd,t,d) is defined in (28). 

Proof: See Appendix D. ■ 






























To sum up, can be evaluated as 



/ qNL \ 

where f ^2 (^) computed by (41) and 

(42), respectively. 


E. The Results ofp™'’ (A, 7 ) and (A, 70 ) 

Considering (23) and bringing together the results from 
Subsections V-A~V-D, p™'' (A, 7 ) for 3GPP Case 1 can be 
evaluated as 

p“'' (A, 7 ) = t \: + (44) 

where Tj', and are computed from numerically 

tractable integral-form expressions using (29), (36) and (43), 
respectively. 

Plugging p“''(A, 7 ) obtained from (44) into (11), we can 
get the result of (A, 70 ) using (10) for 3GPP Case 1. 

Regarding the computational process to obtain (A, 7 ) 
for 3GPP Case 1, only one fold of integral is required for 
the calculation of in (44), compared with three 

folds of integrals for the general case in Theorem 1 . Note that 
an additional fold of integral is needed for the calculation of 
^ASE (A 7 p), making it a 2-fold integral computation. Thus, 
the results for 3GPP Case 1 are much more tractable than those 
for the general case discussed in Section IV, because the linear 
LoS probability function in (20) permits good tractability. 

VI. Simulation and Discussion 

In this section, we use numerical results to establish the 
accuracy of our analysis and further study the performance 
of dense SCNs. According to Tables A.1-3, A.1-4 and A.l- 
7 of [13] and [14], we adopt the following parameters for 
3GPP Case 1: di = 0.3 km, = 2.09, = 3.75, = 

10-10.38, ^NL ^ 10-14.54^ P = 24 dBm, No = -95 dBm 
(including a noise figure of 9 dB at the UE). 


A. Validation of the Analytical Results of (A, 7 ) for 3GPP 
Case 1 

For 3GPP Case 1 studied in Section V, the results of 
p‘^°''(A, 7 ) with 7 = OdB and 7 = 3dB are plotted in 
Fig. 1. For comparison, we have also included analytical 
results assuming a single-slope path loss model that does not 
differentiate FoS and NFoS transmissions [4]. Note that in [4], 
only one path loss exponent is defined and denoted by a, 
the value of which is assumed to be a = = 3.75. As 

can be observed from Fig. 1, our analytical results perfectly 
match the simulation results. Due to the significant accuracy of 
(A, 7 ) and since the results of (A, 70 ) are computed 
based on p“''(A, 7 ), we will only use analytical results of 
p‘^°'' (A, 7 ) in our discussion hereafter. 



Fig. 1. The coverage probability (A, 7 ) vs. the BS density A for 3GPP 
Case 1 with various SINR thi'esholds 7 . 

From Fig. 1, we can observe that the coverage probability 
given by [4] first increases with the BS density because 
more BSs provide better coverage in noise-limited networks. 
Then, when A is large enough, e.g., A > lO^BSs/km^, the 
coverage probability becomes independent of A since the 
network is pushed into the interference-limited region. The 
intuition behind the observation is that with the simplistic 
assumption on the path loss model, the increase in interference 
power is counterbalanced by the increase in signal power in a 
interference-limited network, and thus the coverage probability 
remains the same as A further increases [4]. 

However, the coverage probability performance of the pro¬ 
posed analysis for 3GPP Case 1 incorporating both FoS and 
NFoS transmissions exhibits a significant deviation from that 
of the analysis from [4]. This is because when the distance r 
decreases, or equivalently when the BS density A increases, 
FoS transmission occurs with an increasingly higher probabil¬ 
ity than NFoS transmission. In more detail, 

« When the SCN is sparse and thus noise-limited, e.g., 
A < lOBSs/km^, the coverage probability given by the 
proposed analysis grows as A increases. 

« When the network is dense enough, e.g., lOBSs/km^ < 
A < 10^ BSs/km^, which is the practical range of A for 
the existing 4G networks and the future 5G network [2], 
the coverage probability given by the proposed analysis 
decreases as A increases, due to the transition of a large 
number of interference paths from NFoS to FoS. Hence, 
some interfering BSs may be already so close to the 
typical UE that their signals may start reaching the UE 
via strong FoS paths. 

« When the SCN is extremely dense, e.g., A > 
10^ BSs/km^, the coverage probability decreases at a 
slower pace because both the signal power and the inter¬ 
ference power are FoS dominated and thus statistically 
stable in general. 

It is important to note that the coverage probability perfor¬ 
mance of the proposed analysis for 3GPP Case 1 peaks at a 
certain value Aq. Such crucial point can be readily obtained 
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by setting the partial derivative of (A, 7 ) with regard 
to A to zero, i.e., A* = arg | ~ solution 

to this equation can be numerically found using a standard 
bisection searching [17]. In Fig. 1, the numerical results of 
A* are 19.01 BSs/km^ and 16.52BSs/km^ for 7 = OdB and 
7 = 3 dB, respectively. 

Note that similar trends of the coverage probability vs. the 
BS density A were also observed in [10-12], i.e., the coverage 
probability will initially increase with the increase of A, but 
when A is larger than Aq, the network coverage probability 
will decrease as small cells become denser. 

Considering such trend of coverage probability performance 
and looking at the expression of the ASE in (10), we can 
conclude that the trend of the ASE performance for SCNs 
should be complicated, and it will be investigated in the next 
subsection. 

B. Discussion on the Analytical Results 0 /(A, 70 ) for 
3GPP Case 1 

In this subsection, we investigate the analytical results of 
^ase(a^ 7 o) _ 0,3 ,6 dB based on the analytical 

results of p“''(A, 7 ). Our results of A^®®(A, 7 o) are plotted 
in Eig. 2, comparing with the analytical results from [4] with 
7 o = 0 dB. 



BS density X [BSs/km^] 

Fig. 2. The ASE vs. the BS density A for 3GPP Case 1 with 

various SINR thresholds 70 . 

As can be seen from Eig. 2, the analysis from [4] indicates 
that when the SCN is dense enough, e.g., A > lO^BSs/km^, 
the ASE performance increases linearly with A, which is log¬ 
ically correct from the conclusion that (A, 7 ) is invariable 
with respect to A for a given 7 when A is sufficiently large [4]. 

In contrast, our proposed analysis for 3GPP Case 1 reveals 
a more complicated trend for the ASE performance, which is 
highlighted by the sub-figure in Eig. 2. In more detail, 

• When the SCN is sparse and thus noise-limited, e.g., 
A < AoBSs/km^ (Aq ~ 20BSs/km^ in Eig. 2), the 
ASE quickly increases with A because the network is 
generally noise-limited, and thus adding more small cells 
immensely benefits the ASE. 


• When the network is dense enough, i.e., A S 
[Ao, 10^] BSs/km^, which is the practical range of A 
for the existing 4G networks and the future 5G net¬ 
work [2], the trend of the ASE performance is very 
interesting. Eirst, when A G [Ao,Ai] BSs/km^, where Ai 
is another threshold larger than Aq (Ai w 10^ BSs/km^ 
in Eig. 2), the ASE exhibits a slowing-down in the rate 
of growth (70 = 0 dB) or even a decrease (70 = 3,6 dB) 
due to the fast decrease of the coverage probability at 
A G [Ao, Ai] BSs/km^, as shown in Eig. 1. Second, when 
A > Ai, the ASE will pick up the growth rate since the 
decrease of the coverage probability becomes a minor 
factor compared with the increase of A. 

• When the SCN is extremely dense, e.g., A is larger than 
10^ BSs/km^, the ASE exhibits a nearly linear trajectory 
with regard to A because both the signal power and the 
interference power are now LoS dominated, and thus 
statistically stable as explained before. 

In particular, our finding that the ASE may decrease as 
the BS density increases for the practical range of A for the 
existing 4G networks and the future 5G networks [2] indicates 
the significant impact of the path loss model incorporating both 
NLoS and LoS transmissions on the ASE performance. Such 
impact makes a difference for the ASE of dense SCNs both 
quantitatively and qualitatively compared with analyses with 
simplistic path loss models that does not differentiate LoS and 
NLoS transmissions. As a confirmation, note that in Eig. 1, we 
can observe that increasing the SINR threshold 7 from 0 dB to 
3 dB will accelerate the decrease of the coverage probability at 
A G [Aq,Ai] BSs/km^ because of the more demanding SINR 
requirement, which in turn causes the decrease of the ASE 
(70 = 3,6dB) at that range of A in Eig. 2. Note that our 
conclusion is made from the investigated set of parameters, 
and it is of significant interest to further study the generality 
of this conclusion in other network models and with other 
parameter sets. 

To sum up, our results are different from those in existing 
studies [10-12] and those in [4] assuming a simplistic path loss 
model that does not differentiate LoS and NLoS transmissions. 
The implication is profound, which is summarized in the 
following remarks: 

• Remark 1: Erom the investigated set of parameters, we 
find that when the density of small cells is larger than a 
threshold Aq, the ASE may suffer from a slow growth or 
even a decrease as the BS density increases because of 
the quick decrease of the network coverage probability. It 
is of significant interest to further study the generality of 
this conclusion in other network models and with other 
parameter sets. 

« Remark 2: Our finding is of significant importance for 
practical SCN deployments in the following two aspects: 
fi) The valley area, where the BS density A is tens to 
hundreds of BSs/km^ and the ASE may suffer from a 
slow growth or even a decrease, stands in our way of 
the evolution from 4G to 5G. Since A has been estimated 
as several to tens of BSs/km^ in 4G [13], [18] and tens 
to hundreds or even thousands of BSs/km^ in 5G [2], 















10 


how to cost-efficiently march cross this undesirable valley 
area using new technologies is crucial for the commercial 
success of future 4G/5G networks, (ii) Our results are 
not obtained on conditions such as near-held path loss 
exponents, i.e., < 2 or < 2. All the parameters 

are practical ones recommended by the 3GPP standards 
for the state-of-the-art SCNs. Therefore, our conclusion 
can provide some valuable guidance for operators in their 
quest of network densihcation. 

• Remark 3: The ASE will grow almost linearly as the BS 
density increases above Ai, which might be a candidate 
BS density threshold for characterizing the ultra-dense 
SCNs in future 5G networks [2]. 

C. Discussion on Various Values of for 3GPP Case 1 

As discussed in Subsection VI-B, it is of interest to further 
study the generality of our conclusion in other network models 
and with other parameter sets. In this subsection, we change 
the value of from 2.09 to 1.09 and 3.09 to investigate the 
impact of on our conclusion. In Fig. 3, the analytical results 
of (A,7q) with 7 o = OdB and various values of are 
compared with the results from [4] (a = 3.75). 



Fig. 3. The ASE vs. the BS density A for 3GPP Case 1 with 

SINK threshold 70 = 0 dB and with various LoS path loss exponents . 


As can be seen from Fig. 3, the slow growth or the decrease 
of ASE at A G [Aq, Ai] BSs/km^ (Aq « 20BSs/km^ and Ai Ri 
lO^BSs/km^) is more obvious when the difference between 
the NFoS and the FoS path loss exponents, and a^, is 
larger. This is because the transition of interference from NFoS 
transmission to FoS transmission becomes more drastic. For 
example, when takes a near-held path loss exponent such 
as 1.09 [10], the decrease of ASE at A G [Ao,Ai] BSs/km^ is 
signihcant, and it hardly recovers after Ai. Such observation 
is in agreement with that in [10] when < 2. However, 
the difference between our results and those in [ 10 ] is that 
we have shown how the ASE may also decrease even when 
> 2. Such difference is due to the different dehnition of 
ASE in [10] and in our paper. Note that in [10], a deterministic 
rate based on 70 is assumed for the typical UE, no matter what 


the actual SINR value is. Our dehnition of ASE in (10) is more 
realistic due to the SINR-dependent rate, but it is more difficult 
to analyze, as it requires one more fold of numerical integral 
compared with [10]. To sum up, we draw the following remark 
from the discussion of this subsection. 

• Remark 4: We observe that the behavior of the ASE 
depends on the characteristics of the LoS and the NLoS 
path loss functions. The larger the difference between 
the LoS and the NLoS path loss exponents, the more the 
ASE suffers at A G [Aq, Ai] BSs/km^ due to more drastic 
transition of interference from the NLoS transmission to 
the LoS transmission. 


D. Investigation of 3GPP Case 2 


In this subsection, we investigate the ASE performance for 
3GPP Case 2, which has been discussed in Section V. Note 
that the parameters of the LoS probability function Pr^ (r) for 
3GPP Case 2 are = 0.156 km and i ?2 = 0.03 km [13]. Due 
to the complicated expressions of ( 22 ), numerically tractable 
integral-form expressions for (A, 7) like those in (23) 
are difficult to obtain. Instead, we directly apply numerical 
integration in Theorem 1 to evaluate the analytical results for 
3GPP Case 2. It is important to note that as discussed in 
Section IV and Section V, 3GPP Case 1 only requires one and 
two folds of integrals to obtain p“'' (A, 7) and (A, 70), 
while 3GPP Case 2 requires three and four folds of integrals 
to compute p™'' (A, 7) and A'^®® (A, 70). 

Besides, in order to show the versatility of the studied 3GPP 
Case 1 with the linear LoS probability function shown in (20), 
we propose to use a piece-wise linear LoS probability function 
to approximate the complicated LoS probability function of 
3GPP Case 2 given by (22). Our proposed approximation of 
the LoS probability function for 3GPP Case 2 is defined as a 
3-piece linear function as 


Pr^ (r) 


1 , 


r—di 
d2 — di ’ 


0, 


0 < r < di 
di < r < d2 
r > d2 


(45) 


where di and ^2 are set to 0.0184 km and 0.1171 km, 
respectively. Note that di is chosen as 0.0184 km because 
Pr^ (0.0184) = 0.999 « 1 in (22). Besides, the value of d 2 
is chosen as 0.1171 km because Pr^ (r) in (45) needs to go 
through point , 0.5), which is the point connecting the 
two segments in Pr^ {r) of 3GPP Case 2 given by (22). 


Fig. 4 illustrates Pr^ (r) as defined by 3GPP Case 2 in (22) 
and as approximated by a piece-wise linear function shown 
in (45). For clarity, the combined case with the path loss 
function in (19) and the 3-piece LoS probability function in 
(45) is referred to as the Approximated 3GPP Case 2. Note that 
the accuracy of our approximation can be easily improved by 
fitting the LoS probability function with more than 3 pieces. 
Based on Theorem 1, we can readily extend the results in 
Section V to analyze the Approximated 3GPP Case 2 in a 
tractable manner. The details are omitted for brevity. 
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Fig. 4. The LoS probability Pr^ (r) vs. the distance r for 3GPP Case 2 [13]. 

The results of A^®®(A, 7 o) are plotted in Fig. 5. As can 
be seen from it, the results of the Approximated 3GPP 
Case 2 match those of 3GPP Case 2 well, thus showing the 
accuracy and the usefulness of our analysis with the linear 
LoS probability function in Section V. More importantly, all 
the observations in Subsection VI-B are qualitatively valid for 
Fig. 5 except for some quantitative deviation. 


on the ASE performance of SCNs. Such impact is not only 
quantitative but also qualitative. Previous results predicted that 
the ASE should monotonically grow with the increase of the 
BS density. Instead, our theoretical analysis concludes that 
when the BS density is larger than a threshold, the ASE may 
suffer from a slow growth or even a decrease as the BS density 
further increases. The intuition behind our conclusion is that 
when the density of small cells is larger than a threshold, the 
interference power will increase faster than the signal power 
due to the transition of a large number of interference paths 
from NLoS to LoS. Therefore, our results show that the BS 
density matters regarding the ASE performance, thus shedding 
valuable insights on the deployments of future dense SCNs. 

As our future work, we will consider other factors of 
realistic networks in the theoretical analysis for SCNs, such 
as practical directional antennas. Another future work is the 
introduction of Rician fading or Nakagami fading in our 
analysis because the multi-path fading model is also affected 
by the LoS and NLoS transmissions. 


10 ' 


10 ” 


The proposed analysis [Approx. 3GPP Case 2.,y^=0 dB] (Analytical) 


10 ” 


-The analysis in [4] [a=3.75, y^=0 dB] (Analytical) 



BS density A. [BSs/km ] 


Fig. 5. The ASE ylASE vs. the BS density A for 3GPP Case 2 with 

various SINR thresholds 70 . 


Note that 3GPP Case 2 was treated in [12] by approximating 
Pr^ (r) in (22) as an exponentially decreasing function. How¬ 
ever, the results in [ 12 ] are less tractable than our approxima¬ 
tion. The former results require three and four folds of integrals 
to obtain p^°'' (A, 7 ) and while our approximation 

only require one and two folds of integrals to obtain p^°'' (A, 7 ) 
and (A,7o), thanks to the good tractability provided by 
the linear LoS probability function. 

In the following, we summarize our study in this subsection. 

• Remark 5: The approach of approximating a LoS 
probability function as a piece-wise linear function and 
invoking the results derived in Section V, can deal with 
complicated path loss models in a tractable manner, 
which shows the usefulness of studying 3GPP Case 1 
with the linear LoS probability function. 

VII. Conclusion 

In this paper, we show that a path loss model incorporating 
both LoS and NLoS transmissions has a significant impact 
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Appendix A; Proof of Theorem 1 


Eor clarity, we first summarize our ideas to prove Theo¬ 
rem 1. In order to evaluate p‘^°'' (A, 7 ), the first key step is to 
calculate the distance PDEs for the events that the typical UE is 
associated with a BS with a LoS path or that with a NLoS path, 
so that the integral of Pr [SINR > 7 ] can be performed over the 
distance r. The second key step is to calculate Pr [SINR > 7 ] 
for the LoS and the NLoS cases conditioned on r. 

Prom (7) and ( 8 ), we can derive p“'' (A, 7 ) as 


p“'' (A, 7 ) = f Pr [ SINR > 7 1 r] /fl (r) dr 
J r>0 


I 

I 


Pr 


r>0 

di 


Pr 


'PC{r)h 
_ Ir + Nq 

PC\ jr) h 
Ir + -No 


> 7 


> 7 


/« (t) dr 

/it W dr 


= ^(rL + Tr), 


+ 

+ 

N 


pdi 

/ ^ 

nOO 

J dM^l 

nOO 

J dN-1 


(r) h 

Ir + Nq 

PCh ir) h 


> 7 


Ir + -N( 


> 7 


0 


PCTir)h 

Ir + Nq 


> 7 


fRpir)dr 

/fl.N (’’) dr 

/flTV ir) dr 


(46) 


n=l 


where /]| ^ (r) and (r) are the piece-wise POP of the 
RVs and where and are the distance that 
the UE is connected to a BS with a LoS path and the 
distance that the UE is connected to a BS with a NLoS 
path, respectively. Note that the two events that the typical 
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UE can be successfully covered by a BS with a LoS path 
and that with a NLoS path are disjoint events. Hence, the 
coverage probability is the direct sum of these two proba- 

and 

are piece-wise 


bilities, which should be computed based on ^ (r) 
(^)’ respectively. Moreover, and 


functions defined as " Pr 

and = I 


dn — 


Pr 


PC^{r)h 


/r-t-tVo 


> 7 


PCn{r)h 
R+No 
J-NL 

J 


R n (^) respectively. 


pL 


as 


located at distance R^. 
can be derived as 

(r) = Pr [i?L > r] 

(a 


The CCDF of denoted by (r), 


{Pr[R^>r\X^]} 


r+oo 

I Pr [i?^ > r| = x\ fx {x) dx 

Jo 


(b) 


+ 00 


= / Oxfxix)dx+ Pr [C^^\ X^ = x\fx (x) dx 


>7 /fl,„(r)dr 


^-1-00 

/ Pr [ \X^ = x] fx (x) dx, 

J r 


(50) 


Besides, do and d^r are respectively defined as 0 and oo. 

All the ^ (r) and (r) are stacked into fn (r) shown 
in the step (a) of (46), and fa (r) is formally defined in (47) 
shown on the top of next page, which takes a similar form 
as (3). Since the two events that the typical UE is connected 
to a BS with a LoS path and that with a NLoS path are 
disjoint events, we have (r) = /)j „ (r) + (r), Vn e 


where E[x] {•} in the step (a) of (50) denotes the expectation 
operation taking the expectation over the RV X and the step 
(b) of (50) is valid because (i) when 0 < a; < r, it is 
apparent that Pr [i?^ > r\ X^ = x] =0; and (ii) when x > r, 
the conditional event [i?^ > r| X^ = x\ is equivalent to the 
conditional event [ | X^ = x] . Taking the derivative of 

(l — (?’)) with regard to r, we can get the PDE of Rf as 


{1,..., n} and Y.u=i /d^-i = 1' 

In the following, we show how to compute ^ (r) in (46). 
To that end, we define two events as follows. 

. Event B^: The nearest BS with a LoS path 
to the typical UE, is located at distance X^. 

According to [4], the CCDE of X^ is written as 
x) = exp (—Pr^ (m) 27rMAdu). Taking the 


/L (r) = Pr ^ fL _ 


(51) 


derivative of (l — (x)) with regard to x, we can get 

the PDE of 


fx (x) = exp J Pr^ {u) 2 TTuXdii^ Pr^ (x) 27rxA. 

(48) 

• Event conditioned on the value of X^: Given 
that X^ = X, the typical UE is then associated with 
such BS at distance X^ = x. To make the typical UE 
associated with the LoS BS at distance X^ = x, such 
BS should give the smallest path loss (i.e., the largest 
C (x)) from such BS to the typical UE, i.e., there should 
be no BS with a NLoS path inside the disk centered 
on the UE with a radius of xi < x to outperform 
such LoS BS at distance X^ = x, where xi satisfies 
xi = arg{^^^(xi) =C^(x)}. According to [4], such 

Xi 

conditional probability of on condition of X^ = x 
can be computed by 

Pr [C^^l X^ = x] = exp ^J (l — Pr^ (it)) 2Tiu\dv^. 

(49) 

Note that Event guarantees that the path loss value (r) 
associated with an arbitrary LoS BS is always smaller than that 
associated with the considered LoS BS at distance x. Besides, 
conditioned on X^ = x. Event guarantees that the path 
loss value (r) associated with an arbitrary NLoS BS is 
always smaller than that associated with the considered LoS 
BS at distance x. 

Then, we consider an unconditional event that the typical 
UE is associated with a BS with a LoS path and such BS is 


Considering the distance range of dn-i < r < dn, we can 
extract the segment of ^ (r) from (r) as 

fR,n (r) = exp ^ (1 - Pr^ (u)) 27ruAdu^ 

X exp J Pr^ (u) 2'Ku\d'u^ 
xPr^ (r) 27rrA, (d„_i < r < dn), (52) 

where xi = arg (xi) = (x)}. 

ri 

Having obtained ^ (x), we move on to evaluate 
in (46) as 

7(^r + A'o)' 


Pr 


Ir + No ^ ' 


Pr 


PCk (r) h 


-Nn 


> 7 


= E[7„ 


Pr 


Fh 


h > 


PQi (r) 

1 (Ir+No) 


(53) 


PCkir) 

where Fh (h) denotes the CCDE of RV h. Since we assume 
h to be an exponential RV, we have Fh (h) = exp (—h) and 
thus (53) can be further derived as 

'PCllir)h 


Pr 


Ir 4 

= Er 


No 


> 7 


exp - 


7(/r+iVo) 

P(k(r) 


= exp - 


= exp - 


jNq 

PCk (r) 

JNq 

PCk (r) 


E[/,. 


exp - 


7 


PCkir) 


Ir 


7 


PCk 


(54) 


where .if/,, (s) is the Laplace transform of R evaluated at s. 

Next, we discuss the computation of f^^ (x) in (46). 
Similar to the process to obtain (x), we also define two 
events as follows. 

• Event The nearest BS with a NLoS path to the 
typical UE, is located at distance Similar to (48), 
the PDE of is given by 


i^) = J (l ~ (^)) 2TruXdu^ 

X (l — Pr^ (x)) 27rxA. 


(55) 
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/fl.i (r) 


Jr {r) = < 


Jr, 2 {r) 


I fki^r), 
\ .mr), 
I /fl.2W: 
1 &r), 


when the UE is associated with a LoS BS 
when the UE is associated with a NLoS BS 
when the UE is associated with a LoS BS 
when the UE is associated with a NLoS BS 


0 < r < c?i 
di < r < d 2 


(47) 


fR,N (r) 


when the UE is associated with a LoS BS 
fkkv('^) ’ when the UE is associated with a NLoS BS ’ 


• Event conditioned on the value of Given 

that = X, the typical UE is then associated with 
such BS at distance = x. To make the typical UE 
associated with the NLoS BS at distance = x, such 
BS should give the smallest path loss (i.e., the largest 
((r)) from such BS to the typical UE, i.e., there should 
be no BS with a LoS path inside the disk centered on 
the UE with a radius of X 2 < x, where X 2 satisfies 
X 2 = arg ( 0 : 2 ) = (x)}. Similar to (49), such 

X2 

conditional probability of on condition of = x 
can be written as 

Pr [C^l X^^ = x] = exp ^ J Pr^ (u) 27ruAc?M^. (56) 

Then, we consider an unconditional event that the typical 
UE is associated with a BS with a NLoS path and such BS is 
located at distance Similar to (50), the CCDE of 
denoted by (r), can be derived as 

E-nl ^ Pj. 


p+oo 

L ' 


can get the PDF of as 

rNL 


fk^{r)=Pr[C^\X^^ = r]f^^{x). 


Pr 


Xn 


> 7 


IF] i exp 


= exp - 


pe(r) 

7-^0 \ f 7 




peMo; ^Kpc^^ir), 


Appendix B; Proof of Lemma 2 

Based on the considered UAS, it is straightforward to derive 
(s) in the range of 0 < r < di as 

(s) 

= {exp (—s/r)| 0 < r < di} 

= S exp -s ^ P/3^g^ 0 < r < di 

[ \ ie-Sy/bo / 

*■= exp^—27rA J (l—Ej^jjexp {—sPf3 (u) g)})udu 


0 < r < dij, 

( 61 ) 


Pr [ I X’^L = x] (x) dx. (57) 


Taking the derivative of (l — (x)) with regard to r, we 


where the step (a) of (61) is obtained from [4]. 

Since 0 < r < di, E[g] {exp {—sPf3 (u) g)} in (61) should 
consider interference from both LoS and NLoS paths. Thus, 
(s) can be further derived as 

(s) 

= exp [ —27rA / ( 1 — 


udu 


ro 


1 - E[g] |exp (^-sPA^u I 


(58) 


I y U 

X exp —27rA / — 

\ in di 

X 1 - E[g] {exp } 

X exp 27rA J 1 


Considering the distance range of d„_i < r < dn, we can 
extract the segment of (x) from (x) as 

f^n (r) = exp Pr^ (u) 2 TTuXdu^ 

X exp r (1 — Pr^ (u)) 27r'uAdM^ 

X (1 - PrJ; (x)) 27rxA, {dn-i < r < dn), (59) 

where X 2 = arg (X 2 ) = (x)}. 

r2 

r PC^^(r)h 1 

Similar to (54), Pr i+Mg ^ 7 be computed by 

lilr + No) 


1 - E[g] {exp (^-sPA^^u I 


udu 


.(60) 


= exp —27rA 
X exp 27rA 

X exp { — 27rA J 


1 - 


dij 1 + (sPAL) 


t _ .rwL 


udu 


du 


di 


U 


u 


n it l + (sPANL)”^ 
u 


du 


du 


(62) 


Our proof of Theorem 1 is completed by plugging (52), 
(54), (59) and (60) into (46). 


1 + (SPANL)”^U“''‘' 

to (62), and with s( 
manipulations considering the definition of pi (a, /3, t, d) and 


Plugging s = into (62), and with some mathematical 
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P 2 {a,P,t,d) in (27) and (28), we can obtain 
shown in (26), which concludes our proof. 


Appendix C: Proof of Lemma 3 

Following the recipe employed in Appendix B, we consider 
interference from both LoS and NLoS paths, and we can derive 

7r° 


■^ir I p’anl ) in the range of 0 < r < yi as 




jr 




PANL j 


= exp — 27rA / ( 1 — 


1 - E[g] <'exp I 


iL„ —o 




/ rdi 

X exp — 27rA / 


udu] 


u 


l-E[,]<'exp(-^PA^S-“ 5 


xexp(—27rA / 1 

Jdi 

1-E[g]|exp 

a-i 


udu] 


7r 


PA 




= exp 


—27rA 


v 


1 + 


X exp 


X exp 


/ 

-27rA 

V 

/ 

-27rA 

V 


( „NL 


-1 


udu] 


-du] 


pdi 

I U 

Jr ^1 


l+(irS^^^AlNL)’ 


-du] 


nCO 

J di 


-du , 


(0 < r < j/i). (63) 


For 


in the range of yi < r < di, we consider 
interference only from NLoS paths and we can derive it as 




7r 


^p^NL 

I 

= exp 


X exp 


—27rA 

/ 

V 


Jr dl Y 


( „NL \ —1 

&7PAlNLj . 


-du 


—27rA 


f 

Jdi 


( „NL \ —1 


-du 


iVi <r < di). (64) 


Our proof is thus completed by plugging (27) and (28) into 
(63) and (64). 


Appendix D; Proof of Lemma 4 

Following the recipe employed in Appendix B, we con¬ 
sider interference only from NLoS paths and we can derive 

■^ir ( p7l ) in the range of r > di as 



( NL 

jr 

~~PA^ 




= exp —271A 


nOO 

Jr 1 


-du) , (r > di). (65) 


Our proof is thus completed by plugging (28) into (65). 
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